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MAGNETIC TRACKING SYSTEM 

BACKGROUND OF THE INVENTION 

The present invention relates to magnetic tracking systems of the type wherein a 
magnetic field is established in a relevant work area, and one or more magnetic field 
sensors are operated to sense values of a local magnetic field and are processed to 
determine the position of a tool, instrument or other identified object In general, such 
systems operate using a field generating element or assembly, and a field sensing element 
or assembly, often in conjunction with other position determining elements, such as 
robotic stepper or optical tracker assemblies to track the relative changes in position 
between one or more fixed points or structures in Ihe physical environment or associated 
ixnages, and one or more moving or non-visible points or structures in the work arena. 

Magnetic field generating or sensing assemblies for tracking may be implemented 
in various ways, wilh conventional analog wire coils forming current loops or paths, or 
with semiconductor or microlithographically-formed conductive leads or circuit board 
traces forming current paths, arranged in an appropriate geometry to generate or sense the 
desired field components. TTiere may be a symmetry or duality between the generating or 
sensing elements. Hius. for example in many cases it is possible to have a small multi- 
coil array that generates a spatially distributed magnetic field and a similar or even 
identical array that senses the field so generated. Small coils offer the prospect of 
generating, to a close approximation, dipole fields, although small size may limit the 
attainable field strength or the achievable level of detection signal amplitude. The 
generating and sensing constructions may alternatively employ different scales, for 
example, with relatively large and/or high current coils to establish magnetic field 
components along different axes, and smaller, or more localized coil assemblies for 
sensing field values. Smaller coils, whether for sensing or generating, may. for example, 
be fastened to the body, or attached to workplace or surgical instruments, or to catheters 
or other body-inserted devices, to sense the magnetic field and track position of the 
attached structure.' 



-1- 



In general, it is ttie aim of such magnetic tracking assemblies to define the spatial 
coordinates (e.g., position and orientation coordinates, either absolute or relative) where 
the movable magnetic assembly is located at a given instant in time. It is therefore 
necessary to characterize the magnetic field distribution or signal values with some 
degree of accuracy, and also necessary to accurately detect the field. The field 
distribution may be determined by a combination of field modeling and empirical field 
mapping. The latter, for example, may be carried out as a calibration or an initialization 
step, and may be performed to correct a theoretical field distribution for the presence of 
interfering structures. In any case, the spatial coordinates are generally computed for one 
magnetic assembly (transmitter or sensor) with respect to the other magnetic assembly 
(sensor or transmitter). Typically, one of these assemblies is itself fixed. 

One area in which magnetic tracking has been usefiil is the area of image guided 
surgery. Typical image guided surgical systems acquire a set of images of an operative 
region of a patient' s body, and track a surgical tool or instrument in relation to one or 
more sets of coordinates, e.g., spatial coordinates of the surgical work arena, the 
coordinates of the images themselves, or a target feature of the patient's anatomy. At the 
present time, such systems have been developed or proposed for a number of surgical 
procedures such as brain surgery and arthroscopic procedures on flie knee, wrist, shoulder 
or spine, as well as certain types of angiography, cardiac or other interventional 
radiological procedures and biopsies. Such procedures may also involve pre-operative or 
intraoperative x-ray images being taken to correct the position of and refine the display 
of. or otherwise navigate a tool or instrument involved in the procedure in relation to 
anatomical features of intwest. 

Several areas of surgery have required very precise planning and control. Such 
tracking is useful for the placement of an elongated probe, radiation needle, fastener or 
other article in tissue or bone that is internal or is otherwise positioned so that it is 
difficult to view directly. For brain surgery, stereotactic frames may be used to define the 
entry point, probe angle and probe depth to access a site in the brain. Furthermore, many 
of the foregoing techniques may be used in conjunction with previously-compiled three- 
dimensional diagnostic images such as MRI, PET or CT scan images to provide accurate 
tissue coordinates to which the tracked physical elements may be referenced. Such 
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systems offer valuable improvements for procedures such as placement of pedicle screws 
in tbe spine, where visual and fluoroscopic imaging directions cannot capture the axial 
view that would be required to safely orient the insertion path through bony material 
close to the spinal cord. 

When used with existing CT, PET or MRI image sets, the previously recorded 
diagnostic image sets, by virtue of their controlled scan formation and the spatial 
mathematics of their reconstruction algorithms, define a high precision three dimensional 
rectilinear coordinate system. However, even when provided with such reference images 
it is necessary to correlate and fit available measurement and views and anatomical 
features visible from the surface, with features in the diagnostic images and with the 
external coordinates of the tools being employed. This is often done by providing 
implanted fiducials. by adding external visible or trackable markers that may be imaged, 
and by enabling a surgeon or radiologist to use a keyboard or mouse to identify fiducials 
or features in the various images, and thus map common sets of coordinate registration 
points in the different images. Given a fit of spatial points to image points, software may 
then track changing positions in an automated way (for example, simply transforming the 
coordinates that are output by an external coordinate measurement device, such as a 
suitably programmed off-the-shelf optical tracking assembly.) histead of correlating 
image positions of a set of imageable fiducials appearing in fluoroscopic or CT images, 
such systems can also operate with simple optical tracking, employing an initialization 
protocol wherein the surgeon touches or points at a number of bony prominences or other 
recognizable anatomic features in order to define the external coordinates in relation to 
fte patient anatomy and to initiate software tracking of those features. 

For such applications, electromagnetic tracking offers the advantage that its 
position-defining field, a magnetic field, penetrates the body without attenuation or 
change so that tracking may continue during a surgical procedure, unimpaired by the 
blocking that occurs with visible light trackers (e.g., due to operating room personnel 
moving into positions that obstruct the line-of-sight paths required by optical trackers). 
Optical or ultrasonic tracking, by contrast, may require a larger or excess number of line- 
of-sight paths, and corresponding transponders and/or detectors to assure that 
triangulation is always possible despite occluded pathways. The body-penetrating 
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electromagnetic fields also allow one to track locations or movements inside the body 
with minimal resort to the fluoroscopic or ultrasound techniques normally required for 
visualization. ' 

Among electromagnetic tracking techniques, several principal approaches are 
known In one, relatively large Hehnholtz coils establish well defined and highly 
uniform independent magnetic field components or gradients along each of the ^. 7 and Z 
axes in a work arena, and the static field components are detected by a localized detector 
to determine position coordinates. Ihis approach has been proposed, for example, for 
cranial surgery, where such coils may define a suitably localized region encompassmg the 
entire operative arena. Another principal approach involves using time-varying dipole 
fields e g., dipole fields establishedby driving field-generating coils with an AC current 
signal While the latter approach offers some processing advantages (such as being able 
to synchronously demodulate induced signals, and thus cumulate detected signal values to 
enhance sensitivity, and also the ability to establish &e^. 7 and Z field components at 
different frequencies so that detected sensor output signals may be separated or even 
demodulated simultaneously), it has the disadvantage that varying magnetic fields mduce 
eddy currents in conductive structures found within the field. Induced currents 
themselves then generate secondary magnetic fields, thus introducing distortions into the 
expected or calibrated field distribution. Conductive or metal stmctures are in fact 
commonly present in a tracking enviromnent. whether it be an avionics, surgical or 
industrial tracking environment 

The latter problem has historically been addressed by the observation that for a 
fixed metal disturber located at least twice as far firom the field generator as is the field 
sensor the induced disturbance will be low. e.g.. under one percent, so that by restricting 
the tracked arena to a region sufficiently removed from the disturbing structure, accuracy 
may be achieved. However, this approach may be unrealistically restrictive for certam 
applications, including some image guided surgery applications, where highly disturbmg 
equipment (such as the imaging assembly of a fluoroscope) is necessarily placed as close 
as possible to the work arena in which tracking is to occur. Another common alternative 
approach would be to map the disturbed fields or detected signal levels that occur close to 
the distorters present in the work arena, so that a processor can more accurately determme 
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^rdinates torn fte r»„ toe field values or ftc induced signrfs de«c«d by . sensor. 
However, as a practical matter, such mapping is no. only likely «> require . ttne- 
consumingpreliminaryse.-upoperation.bu.may require thar^eset-upbecarrredou. 

aftesh for each new arrang=men.ofoperating room equipment .ha. in«odueesdtfer«,. 
distortions, i.e., when equipment is changed or moved. 

to addition u> the practical problem of accurately detecting to field, there is .he 
fceoretical probl«. of converting the signal measurement, into position and oricnta«on 

coordinates. 

On a fundamentd level, fte task that must be addressed by any electromagnetic 
«ckingsystomistocompuU«oual.yd..ennineauni,uepositionftom*eva™^^ 
measured parametotsCWi-nyinduc^lvohages indicative of field stiengtir). Often fte 
„,ev».. equations have several solutions, and care mus. be Uken to operate w,ft,n a 

sir^e-valued domain (Wi-lly a hemisphere, quadran. or oCan.). ftus limitog fte 
seitionoffteimtia. generator orreceiverfixedlocationstoes.ab.ishas„ffic,entiy well 

behaved field region. ^ restticting fte allowable work arena in relation ftereto. or else 

providing «lditional, or «c.rinsic data inputs to resolve ambiguities or refine 

companions. Beyond encountering multi-vataed coordinate solutions, fte accuracy of 

ft. cooniinate determination processing may depend quite critically on fte relative 
positions of generating «.d sensing element. Processing «r-«ons may break down or 

solutions ma, become poorly defined as fte sensing or generating elemen. approaches or 

crosses a particular plane or axis, or one of iB piteh. yaw or roll angle coordmates 
approaches 0 or 90 degrees, becomes too acute, or becomes too obhque. 

Anofter practical problem stems fiom fte dynamic range of fte sensed signals, 
which varies greatly wift disbnce^sition, and may drop to quite low levels wrth 
increasing distence from fte ««.smitter, or wift decreasmg size of fte genemting or 
sensing coils. Often, as a field mn. moves to differen. positions, different gains mus. be 
determined on fte fly. or additional gain steges mus. be used in order to obtiun adequa.e 
signal values. This intioduces some complexly- and potential for error in nom«Imng fte 
signals to accurately fi. .ogefter readings fiom two different or even closely contiguous 
regions. 
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Thus, while magnetic tracking offers certain significant advantages, particularly 
for surgical applications, there remains a need for improved systems and processmg to 
enhance accuracy of coordinate determinations. 

Accordingly, it would be desirable to provide an electromagnetic tracker of 
enhanced accuracy. 

It would also be desirable to provide an electromagnetic tracker having enhanced 
immunity to common field distortions. 
SUMMARY OF THE INVENTION 

One or more of these and other desirable features are obtained in an 
elec^magnetic tracking system wherein a field generating unit and a field sensing umt 
are arranged to generate and to sense, respectively an electromagnetic field in an aren^ 

interest. At least one of the units is movable, and the units are com^ected to signal 
conditioning and processing circuitry that detects the levels of*e transmitter to^^ 

signals and ihe received signals, ratiometrically combining them to form a matrix 
representative of the mutual inductance of each office pairs of component coils. The 
n^utual inductance information, providing functions of flie relative positions and 
orientations of Ihe two units, is then processed to determine corresponding coordmates. 

Since mutual inductance is a symmetric property, i.e., depending simply on the 
geometries of flie transmission and receiving coils and their relative positions, the system 
xnay interchange transmitter and receiver units, employing quite small coil assembhes. 
yet useamodel that givesadirectcomputationofcoordinates without excessiveiterative 

approximations, and without resort to the multi^de of gain level and other normalization 
corrections otherwise typically needed when working from magnetic field intensity 

measurements. The equations may be processed and solved in real time to effect surgical 
or other position tracking, and various corrections and calibration of signal magnitude 
(e.g.. for cross-coupling of coils or for finite-size and non-circular coil geometries) are 
readily implemented. 



to accordance witt, anoto aspect of *e invcdon, 4« m.g=«ie nactog sys«ns 
of fte invention may employ signal conditioning or processing electtonics m winch 
xelev^.. signals pass toough fixed gain amplifies u. high precision (e.g, 24 ba WD 
digiW converters. This bi, size encompasses a dynamic range effective for mul«-b.. 
^Ution of signals enco„n«ed a. all regions of *e intended «ackmg «ea, «,e^ 
aLing *. need for variable gain or AGC element, or for multiple or drfferen, preamp. 
having difibren. gain levels, as one coil assembly moves firrdrer ftom fte ofcer. Thus. *e 
provision of a common preamplifier with a high precision sampler eliminates patehmg 
and fte errors due to inaccmate gain ratios between gain states. 

ta accordance with another or ftrrther a^ of the tov«flion. the processor 
employs a model or se, of equations *at open,te. with the sensor output values te 
determine the position/orientation coordinates of the «nsmit«.^ver assembhes, ^d 
„hen a tracked magnetic coil assembly element approaches a singular region. ..e a plane 
or region where the model becomes ill-defined, inaccurate or laclcs a solution m fte 
coordinate system, the processor operates by tiansfbrming to a coordinate „te«on 
inwhichthedetectedvalueslieinawell-definedornon-singularregion. Theprocessor 
ften solves to determine a coordinate measurement, and transforms back to fte 
underlying or original coordinate system. 

accordance with another asp«:.. an open.ting room system of fte invention may 
employ a conductive shield, or a shea* stiuchrre confi^ to fit about or contain an 
interfering component or piece of equipment The sheath standardizes the magnettc fie d 
disturbance inttoduced by the component In some inst««es the sheath may be a „eUl 
cylinder,dimensioned.oenclosethedisnHbmgpiecofe,uipment Thesheafhor 

cylinder may be formed of sheet metal of a gauge such that eddy currents are mduced by 
the magnetic field, thus giving rise toastandaniized field distirAanceorigmamrgata 

contour or surface external to Ae equipment Altem-ively. ofter suitable condtrctive 
materials, such as a catton fiber composite material, may be used. The sheath also 
moderates or effectively nulls electtomagnetic disterbances originating within its contour, 
eg attheunderlyingpieceofequipmentdisposedwitirinthesheath. The incluston of 
this standardized dismrbance in a system of the invention permits the interchange o 
different pieces of equipment without introducing excessive variations in resultmg local 



disturbance ch.«c«ris«cs. Dis«,rbanc.s r^y tos be mapped, or even n,ode.ed in a 
dngle initial se, up »pe.a«on. without necessitating compiUtion of a new d^o. 
map for each new piece of equipment 

in accordance with another or additional a^t of the invention, mther than 
simply introducing a sheath to fonn a standardized disturbance. ^ ^^^^ 
such a disb^bance. For example, the processor may model a conducbv. sheath fitted 
about a certain re^on as a conductive ring or cylinder a, *.t region (using the tatown 
din,«,sions and behavior characteristics ofthesheetmetal material). Tl»modeled 

distutbance may then be added to the stored values of a map of the undis»rbed 
electromagnetic field to form an enhanced field map. or may be o«,«w,s. appl»d to 
enhance accu^cy of tracVdng detenninations. Tbe modeled field may also be used «. 
provide a seed value for determining position and orientation coordinates. A fitting 
procedure then refines the initial value to enhance the accuracy of ti« P&O 
determination. 

BRIEF DESCRIPTION OF THE DRAWINGS 

-n.ese and other features of the invention will be understood from the description 
below and accompanying claims, taken together with the Figures iUustratmg 
embodiments thereof wherein: 

Figure 1 illustrates basic hardware elements of a tracking system in accordance 

with the present invention; 

Figure 2 is a hardware block diagram of a preferred embodiment; 

Figure 3 is a high level flow chart illustrating operation of a system accordmg to 

one aspect of the invention; 

Figure 4 illustrates signal processing for determining coordinates w,th the system 

are 3; 

Figures 5A-5F illustrate circuit details of some exemplary embodmients; 
Figure 6 illustrates aspects of distortion correction in accordance with other 

aspects of the invention; 

Figure 7 is a network model of the system of Figure 6; 
Figure 8 illustrates geometry of fields for the configuration of Figure 7; 



of Figure 3; 
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DETAILED DESCRIPTION 

The present invention is an improved electromagnetic tracking system which may, 
by way of example, be employed in an operating room or clinical setting to determine the 
position and orientation of an object such as a catheter, surgical tool, marker or the like in 
relation to other physical structures, or in relation to images that may, for example, have 
coordinates assigned thereto by other means, such as via coordinatized 3-D inwge sets, 
position-defining frames, or correlation with imaged or tracked fiducials or markers. A 
great number of such systems are known, being described for example in United States 
Patents No. 5,967, 980 and related patents of M. Ferre et al. No. 5,307,072, and others, so 
the features of such systems will not be further described here. 

As relevant to a first aspect of the present invention, an electromagnetic tracking 
system includes a magnetic field generator and a magnetic field sensor together with 
signal conditioning circuitry and a processor for processing Ihe conditioned signals. The 
generator and sensor are also referred to as transmitter and receiver below. In general, a 
varying current injected into one cofl assembly produces a magnetic field, and this field 
induces a voltage in the other coil assembly. The induced voltage V is proportional to tiie 
rate of change of the applied current I, with a constant of proportionality equal to L„y the 
mutual inductance. According to one aspect of the invention, a magnetic tracking system 
is configured to apply the transmitter drive signals and the induced receiver assembly 
signals to determine mutual inductances, from which the processor flien determines the 
coordinates - i.e.. the position and orientation coordinates, or simply "P&O"- of one coil 
assembly of Ae system with respect to the other coil assembly. 

For two closed circuits the mutual inductance is given by Neumann's formula, 
and is equal to an expression which depends only on the geometry of the closed circuits 
(e.g., the coils), and not on the electronics used to measure the mutual inductance. Since 
fte mutual inductance L„ is a ratio independent of applied current, wave form or 
frequency, it remains a well defined measurement parameter under a wide variety of 
conditions, offering improved performance in real operating environments. With three 



transmitting co«s and ihree receiving coilsTn^oe" mutual inductances preside an ov«;- 
determined system for derivation of the position and orientation coordinates. 

In some embodiments, the roles of transmitter and receiver assemblies may be 
interchanged, each comprising a substantially identical three axis dipole coil assembly. 
Preferably these assemblies are built so the three coils of an assembly each exhibit the 
same effective area, are oriented orthogonally to one another, and are centered at the 
same point. For many of the applications alluded to above, the coils are quite small and 
ideally are small enough (compared to the transmitter-receiver distance) that they may. as 
an initial approximation, be considered to exhibit dipole behavior. However, more 
generally, the system includes embodiments discussed more fully below, wherem one or 
more of the coils are modeled or feeir signals processed as finite size transmitter or 
receiver. 

It should be noted that the closed magnetic generating or receiving circuit includes 
the coil, its comiecting cable and some electronics to which the cable is comiected. In 
accordance with a preferred embodiment of the present invention discussed further 
below, contributions of these extrinsic or ancillary circuit elements are reduced or 
eliminated, permitting a simple and highly accurate model of the coil alone with a direct 
short across the coil leads. 

Figure 1 schematically illustrates a basic measurement system 1 employed in a 
prototype embodiment A transmitter assembly 1 0 is driven by a transmitter driver 1 1 to ■ 
produce a magnetic field, while a receiver assembly 1 2 responds to the magnetic field by 
producing induced signals, coupling the signals induced therein to a receiver voltage 
preamplS and an analog to digital converter 14. Tbe transmitter assembly 10 and 
receiver assembly 12 may, for example, each include three coils, oriented along 
respective spatial axes. The digitized values of the receiver voltages from the digital 
converter 14 pass to a digital signal processing system 1 8 that may be implemented m or 
may interface with a host processor 20 Values of the transmitter coil drive signals 14 are 
also digitized by an analog to digital converter, and pass on line 15. to the digital signal 
processing system 1 8. In a preferred embodiment, the driver 1 1 applies known or 
approximately known voltages, rather than known or approximately known currents, to 
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«,e ™i«er a^ly .0. ^ — »„c n»re induct » -^"f 

tt,e ^ansmitt^^bly. 11,13 votege^v^ approach advantageous^, a..ev,a«. it. 
„.ea «> co„.ol *.e phase of o« or more applied .o *e «nsnu«« co.U and 

6.rt,erpe»teftesys«m«.toe«onevenwid.n»s»neda„d/.,no„^rt.o^nd 

coUs. Those having ordinary s«U in *e ar. wiU apprec a« fl. *e — 
^ a sysun. according .o *,e cachings of *e invendon e«. also he drwen by currents 
rather than voltages. 

various signals are tt.nsmi«ed, sampled, and combined and processed in 
^tfon, for .hich operation a timing and synchronization con». unit 30 fonns a 
numberofbasicsynchro.izedtimingsignals«>deSneorco„«l.forexampi^*e 

mu,tip.exingofsignalsbe.wee.«.~siona.^ receiving coils, conrro^^^^ 

Ae digitizing circaiK conBol analog.U,^gi»l (ADC) samplmg ra«. and fl,e 
processingoftiredigitizedraw signals ftomftecoils^produceaprccessedrawdau 

„u,putThehos„»cessora,s„opera«s«id.*era«da.a«,efreC*ecompaU«onal 
algoritos for determming coordinates and correcting position and onen«tio„ 

determinations. 

Figure 2 shows amore deoiled hardware block diagram of Ms aspect of a 
^g system which, in acconlance with the present invention, is confl^ to m^e a 
robustandaccuratedeterminationofmumalinductanceofthevarioustransmrt^and 

liver assembly coil pairs, m. information is then applied as fi.r.her descnbed beiow 
«> detennine theposition and orien.ationcoordina.es of one assemMy --^^^^l 
other. Forsimplici.y.ti,esame„m.er„sasemp.oyedinFigurel areusedtodes^te 

the transmitter, receiver, host processor, digital signal processor and synchronrzation or 
control circuitry. 

AS shown in this embodiment, a Cock generator circuit 30 provides timing sipials 
which are applied to different units for effecting synchronized processing descm«d 
further below. These units mOude a signal sampling analog-to-digiul converter 1.4. 
(sometimes referred » as ADC 1 14 or simply an A/D unit below), which d.g,tizes ^ 
signal values received from fte receiver assembly and fte drive signal values «>mpled 
fil fl.e —er assembly; a digiral-te-analog converter (DAC, 1 14b whtch converts 
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digtol drive Worn, signals .o analog votog« for powering ft. amplifier »Weh 
Jpplies ^nsmitter drive signals; ».d a digi.al signal p^cessing uni, 1 8 and wave fonn 

generator 40. 

TT,e wave fonn generate, 40 is a digiul ftequency synftesizer to. generates 
digiUl values of a sine wave signal a. one or n,o,e selected frequencies. Ttese n,ay 
inld^ for example, a separate ftequency for driving each of *e transnutte, cotls. mse 
digiul values are converted to analog values by DAC 1 Hb. and are fed to a power 
a^upHfier and filter 50 fi.arivi„g-h«ansnn«er coil. Advantageously.thepow 

antpUfiern^ybeavotageantplifier. Tl. drive signalat any given ..W,,ssw*^^^^^ 
e^L. multiplexer to the appropriate coil, and simultaneous^, the n,duced voluges m the 
^ving coils are demodulated a. the ft«,uency of die transmitted drive s.gnal. 
Meanwhile, thedrivecunentto retransmitter assembly isdeteetedbyacurrentsensor 

5U and is measured by current measuring circuitry 60. 

AS ftHther shown in Figure 2. ^ receiver coil assembly includes a calibration 
t^ttingconiOafixedatapositionclcseto the tee receiver coilsofthc^semb^ 
CoillOam.ybe.fbrexample. innnovably secured in a fixed orientaton rela«ve to the 
teee receiver coils in an integrated unit TTe coiU Oa is operated generate a local 
magnetic field, oriented so that i. induces calibtation voltage signals in each of the 
^vingcoas.Tl,edrivecn,rent in this coil is also mea^tedby the cnrrentmeasunng 

circuifyeo. Coil 10a is energized, as described ftuther below in connection w,th 

detennination of muntal inductances, to effect normaH^tion of .he component stg^ls 

^ed in the tee receiving coiUR..R.R..-n»c.litaati» 

calibration coil driver. 

The receiver coils R„ R,. R, provide to votoge outpuB to a preamp and fB^ 
80 the output of which is fed to the analog-to^igital convener ADC IMa. as described 
ab^ve. to produce digital signal values for processing. Converter 1 14a also recedes as an 
i.,pu, tite outputs of d. cutren. measuring circuitry 60. so to. the ADC U 4a provtdes 
digital values of both to calibration coil current and to transmitter coil currents to to 
digitalsignalprocessingunitl8.Coordinationbemeentovarioustiansm,nmg^ 

receiving, sensing and converting steps is carried out with timing signals from to clock 
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,^,u>, « 30. By way of exan^le. Cook ciro«,«y 30 n«y '^^^^ ^ 
Lgahe«zolocfcdividmgi.dow„(e.g..«>y256).o».a„AIX:sa„,p.«ra«of,7.65 25 

Which is s™ly diviaea .y 2, ^ 1440. and by 256 . p^vde s.n^.0 
Ling igna,sof4S m.25 H,oh«1.for*.DAC.a.d«,dofcyc,efta,n= s.^ a. 

^^^2 ^ for a .™ of produo. (SOP, pr^S n,odu,« i^.en,en«d ^ d» s.g^ 
!;IessorDSP18.Fu,*erfte,u=„cy division ™ypro,id.ao.™..n,ump.«cr cycle 

::Zs„i.hing«ns„i«e,d.veoyc,=s.odiffc™.coils)ofO.,3245477Hc«z.c. 

Jsynctaoni^.in^.si^>--PP-P^- ^us.op.««o„of«,.va„o^ 
Ld,I«ni« Shown inFigu,e2issynchroni.edby«,ec.ockgcnc«.rc™:u,.30. such 

an Komplary timing unit is shown in Figure 5E. 

Oneof*oc,ocksign.lsis.W»adigi«wav.fo™g«.»tor40U».syn*^^ 
a«edis«nc.d.vo^ucncies^*..-,.and..»nsn,i««coiU.«w.n^.ood 

d,a.d,edigi.lwav=fo™.g.ne«».40™yhein,p.cn»n«d,n*ehos.proc^20. 

However.i.isp.efcrabIyaspecia.i.»dft«r-cy^«-»*'P»'^*-7'f' ^ ^ 
plce;adigLcon«nuonssi„=wa,eofpr.ise6«^e«:y--P'-As^-^ 

Ue.below.d,edifreren.fte,nenciesareuscd«> drive indcpendcnt-gnehc field 
rene..,coils,andco„espondingsigna>sinducedinreceive.coilsarede«c^»^ 

p^ i. a sin,p.e digiul demod„.a«.n a sun.-of.produc« processor. *e 
syn*»sized,^e,„enciesarein,egrai™,.iplesof*esnn,ofproduc«d,g..a,s,gnal 

:Lssingcyc,e,a.esCi.e.,*ecyclera«sof*eI>SPnniUS)so,...*e^e^^^^^ 
repea. exacUy every SOP cycle wiO, no discon«nui^ be»-een cycles The ft^»c.es 
so generated are fixed, and «.e sine wave generator n«y be in,pl«nen«d .n a ^ 
Jion by precalculating samples of desired sign wave signals and stonng *e 
lasmonoypi o n5^^mv^^ The PROM addresses are tiien 

c,mr,l^<5 in a OToerammable read only memory (PROM). lUe i-K^^m a 
:rys:l"n..in.inggeneratingconn.ers(no.shown)so.hat.ePK^^ 

sanjles in the s^al fonna. required by the digital-to-analog converter 1 14b, Another 
PROM output may implement the divide by 1440 contwl logic. 

The digiul output values of generator 40 are converted to a differential analog 
voltage in d,eDACcircui.^ll4b. which may.fbr example. bea24.bi.l^«as,gn.^^^^^ 
digit^l-analog-converter. Its output analog signal is amplifiedby the amphfier 50 wh,ch 
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"n«y.forexample.betap.emenWwi*alowpassffl.er3ndafix.dg.indiff«^^^^ 

amplifier. 

AS is appax^,. fiom fhe circuity d«cribed tas far. systems of *e mveutfon sc. 
„„,»„easurebo*U>ettansmi«erooadrivecu™>.sa„d,he«ceivercoiHnd„.ed 

voltages. TT-escaretatiome-ically combined «>fonnmutua.mduc.an..—men«^^ 
Vanl construc^on details ™y be practiced in otder to have me detected st^aU 

as accurately as po^ible the desired positional field dep«.dence. F.garcs 5A- 
5D illustrate fcrfcer construction det«ls useful in preferred con^cttons of the 
transmitting and receiving hardware. 

AsshowninFigure5A.eachtransnu«i„gcoilisdrivenbyadifferential.n,pMer 

50 that antplifics .he analog-converted syn«.esized drive iie,uency signal r of 
the differential antpMer is connected to *e low end of fte ,ransm,«er co, (denoted i m 
the Figure), which is also connected u,avirt„al^ound.Acoaxialcable55cames the 

arive signal to the Ugh end of the transmifler coil I. fl.us avoiding int,oduc«on of 
elecnnnagnedofieldsduetocurrentswiftinthedrivecable. The virtual ground, 
although no. exacfly a ground voltage, is connected to d,e , input of d.e amphfier. *»s 
that the acuta, voltage applied across dte — er coil will no. be affec^d by 
s«ay voteges on fte virtual ground. In this n«nner. a stable drive voltage ts 
applied across Ae transmitter coil. 

Further, because the impedance of amplifier r input is high compared U, to 
impedance of virural ground or real ground, and its current signal is low. fte current 
d^wn through *e r eomtection is substantialV less than (e.g.. typically ,0 mnes less 
to) *. current ftrough the coil it^lf Ms small error may be estima»d and correc«d 
mathematically in subsequent processing. Similarly. *e coaxial cable capact^nce may 
*un,asmallcurren,aro»ndtocoil;d.iscurren.mayalsobees*n.»dandcorrec«d 

mathematically. 

AS further shown in Figure 5 A, as a coil is driven by .he amplifier 50. *e coil 
curren. is sampled between flre low em. of *e coil and the r tenninal of fte dnvmg 
amplifier. The de.ec.ed drive cunen. is applied .ogefter wi* induced recetver co.l 
signals for ttie compuation of muttal inducance. 
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I„ addidon. shown « Figu. SB. a referee in,p«U„c '^"^''^'^ 

rlfer^ce impedance Z.. T.ei.^anceofZ^ismuch.a^er««.meooU.ocoU 

T H so the Z ^ driver provides a smaller current As seen m Frgare 5B, a 
mumal nnpedances. so fl.e Z^ dr. p ^^^^ 

resistor J!=whjoh is large compared to Z,tf. converts roe o s 

Zm«^-.o a small current in Z.^. Ws connection. W... the difference amplrfl^ 

:;^gro„ndea.ra*erthancor.ected.otHelowe,^ofZ,.for.»oreaso^ 

Z inp« r draws a sm^l curren. which U significant compared - " ^ 

tZ. *e receiver measurement channel responses to ^ are me.^ on^ when the 
c.rren.flr.oughZ^isbeingmeasuredbythecurren.-meas«remen.cha«,e..T^^ 

reference in^ce may plug directly onto the connectors on the recetver 
^lmicsWa.us.voiding*ein«duc.onorar.ifac,sorthep.^«ce^cter,s.cs 

, J- ♦ 7 ;c ,,tili7ed only for determining the mutual 
of cabline. In fliis embodiment, Zr^ is utilizeQ omy xui , ^ , 

^:I^LL.h«ween.heca,.-hra«on«nsmi.tinscoillOa(riOUBB.,a^^^^^ 

Livereoils atftedmeofmanuftc^Hng^ereceiver assembly. Wta^^c,.ng..he 
^U- valuesco.responaing.om„.u..i«ductancesbe.ween,^ec^.bran^ 
«^„ing coil 10. a»> each of Are three receiver coils are employed as m„mal 



inductance standards. 



The overall arrangement of the transmitter seeBon is shown in Figure 5C. The 
ariver/amplifler 50 is comtected ta nun to each of the coils and referenced hy a c^.- 
Lh^mu,.iplexerCC^DO60. The virtaal grounds for aH the .ransmrtter eorb go 
nracurrJl-n^deswitehma^thatconnectsnoneoronetransmi.^ — 

to theVrimary winding (e.g.. a single-urm primary) of a current transforrner r„ whrle 
oonnecdngallthe remaining transmitter coilstothehigh-currenttransmrner^.. 

common-ground pota. The selected current .us goes tluough -'J^J ''^ ^ 
current «nsformer n to the ground point A suitable current transforme^ r.n«y have a 
^^-to-secondary ^ ra«o of about 1:200. It operates to isolate d.e U.^—«« 
Lents tom «.a.og ground of the receiver elect^nics. accurately reducmg the 
«nsmittercunen,stovataeswimin.herangeof.hevir^al-groundpreamp.and. 

minonizingthereflec.edvir»al.ground impedance seenbymetransmtttercorl. 

Tic munral-impedance reference is preferably operated a. a lower current than the 
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coils are operated, so 4,e ™^.nnp«l.nocref«e»ce cu™.. is s«.«h«i -» 
or ou. (as shown) on fte secondary side of current transformer T,. Moreover, toe 
secondary of is always in fte circui, u> cancel effects due «> Bansformer magnettzmg 
cu^,. so ,ha. the magnetizing current then appears as jus. another con«an. con,plex 
gain error in *e current-measurement channel. Such constant gdn errors cancel out m 
data-reduction calculations. In Figure 5C. the op amp and resistor «,fonn a v«uaU 
ground current-mpul, voltage-oulput amplifier. 

HGURE 5D provides an arrangement alternative to that shown in FIC3URE SC. 
to to altemative arrangement three transformers (62, 64. 66) are utilized, and fte 
current switches are connected to the secondary windings of the transformers. Tie 
i,„prfa„ce vacation seen by each transformer coil in the transformer primary crcu... as a 
resultofopeningandclosing the corresponding switch, ismuch,owerthan,h.var,ation 

seen in the «rangement of HOTRB 5C. In particular, the variations are reduced hy a 
Uctor equal to tl« scpare of ft. transformer turns ratio. For example, the use of . tun« 
ratio of 1200 in the ammgement of FIGURES 5D results in reducing the tmpedance 
variations by a fector of 40.000 .eUtive to ftat in the arrangement of HGURE SC. 

In operation of system, magnetic fields generated by the transmitter coils 
induce voltage signals in the receiver coils, and these are received and sampled to 
develop raw received signal values fcr the processing. Tbe voltage-inpu, preamp and 
ADC channel fer one receiver coil i is shown in Figure 5E. "The voltage from the 
inductor L is appHed via a coaAl cable 55 ,0 the differential-voltage-input preamp. Tlte 
™mp output is applied «> the input of an analog to digital converter, such as the above 
describ^l 24-bitD.lt.si^ADCthrough.passiv.lowpassfil,er. T^e filter ehmtnaes 
alias responses as well as eflects due to nonlinear ADC input currents. Use of a coaxtal 
cable has the advantageous benefit ftat extemal electromagnetic fields can induce no 
potentials beti.ee. the center conductor and the shield of the coaxial cable (assummg a 
perfectiy-conducting shield). By contrast, a conventional twisted-pair cable may cancel 
effects due to mnform appli«i electromagnetic fields, but can introduce extraneous 
signals in *e presence of gradient fields (such as spatially varying fields) are used m tire 
system. 
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Thus, in *e illus»a«d embodtacn,. tt>« is one voltage input ehanne. for each 
receiver coil, and a channel for the current preamp. 

overall operation 200 of the position sensing system proceeds as shown in Figure 
3 TT„ transmitting.-.^- and z- coils T„ T. are driven to gene^te magnetic fields, 
and the voltages induced in the three sensor receiver coils denoted R., ^ "° 
:Lea.:naddition.areferencemagnetic«eldisgenera.edina — coU.^^ 

eoU toedinoro,osetofl,er«=eivingcoi,assen,b,y. The currents flowmg though -he 
HsmittingcoiU and a,e current flowing ^roughthecalihrationcoa are eachmeasured. 

and an the signals, including the received voluge signals in each receiver cotl. are 
:^.i.edand^sedtoasignalprocessingunitdescrihedhelowmore*i.nymeonn«.^ 

w,^Figure4. -m ADC sampling is controlled by a common timing 
«,e tini;^ relationship between ADC samples and DAC samples is tightly cont^Ue^^ A 
suitable timing module is shown in Figure 5F. having the representanve ftammg, c^de 
::e and o Jp-oeessing or sampHng intervals described above. The ADC outputs, a^ 

^ reference information .iom the timing generator pass to a digital stgna^ proe^smg 
^ signal processing applies ^dn.mre demodulation (i.e.. dtgnal heterc^y™ 

!.hed,iving^encyofthec„rrent.y-ac.ua,ed^smitting.-.^..-orcahb«^^^^^ 

signal) to accumuUte smn-of-products (SOP) complex raw si^l values, and computes 
vSous matrices which allow one to nonnalize the coupling between the cotls and the 
response of the leceiver coils to toe caUbration coU. 

These are used in toe next stitge f form a nine element mumal inductance matitc 
for all ofthe (3) transmitter and (3).eceiver coil pairs. The mafl» is then used to 
compu.ecoordina.esofthereceivercon.ssemblywi,h.especttotoe— er.The 

overallapproachfortoelattcrcomputationmaybesimilartothatemployedma 
conventiL field mode, dipole approach, such as' toat described in the 1,79 paper of 
F H Raab e, oi Volmne AES-,5.Ko. 5.pp. 709-718. andFrederick H. Raab^Quas.- 
static Ma^etic Field Technique for De,«mining Position and Orientation. IEEE 

« • i^rj 1Q r^^- 7'^5-243 October 1981. 
2Ta«J»cftonso-.(3«.science<.»rfJic».o<eSensmg.GE-19(4).235 243.u 

For computation of position and orientation coortina.es. toe processing preferably 
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p^eeds by surting wifl. a. initial estimate of d.e six position and onentaOon 
L.ina«s.app,.i..«.osecoo^ina.es.con.„»an„ti.a.ina.c«n^^^^^ 

i^tively refining fl,e coordinate estimate to enhance the goodness of fit Ween .he 
™.n.al inductance values as measured, and those predicted from .he estimate, 

AS discussed further below, various fe^tirres of the processing make the mutual 
i„ductance model highly accurate and ,»t prone ti> errors of scale or positi^on. Mor^v« 
direct solutions are possible, the iteration proceeds <,uic«y so that the ,^n ^d 
citation coordinates ma, be output based on *e estimated coordinates and me.^ 
lunral inductance matrix. Appendix A attached hereto ilhrstrates a su..ab.e ca^ul ti„u 
I deriving an initial set of position and orientation coordinates fiom a measu^ set of 
dipole-modeledmumalindu«a««L.ofthe.ransmi««:/recewercodp«rs. Th 
lensionsof^coi,sareassumed«>besmal.comparedto,heirs.paration..llowmg 

dra:lrorexample.asdescr.W..hepaperofPredeHc..Ka.^^^^^^ 

on Geos^ence and Remote Sensing. OE.,9 (4): 2^5-243. October 198 . M^i 
Luctances may be calculated for dipole coils or for circuits compnsed of strargh. Ime 
mductancesmy Aooendli C Each of the foregoing appendices 

segments as outlined in Appendix B and Appendm-. 

isherebyincorporatedhereinby reference andmadeapartofthtsdescrrption. 

Figure 4 is a block diagram showing the digital sig«a processing 400 undertaken 
with the dTgitized coil drive currr«,t «k. sensing voltage measurements. As shown, each 
:fr:g„!.measuremen.sisdigitizedinahighprecision(24bi.)ABC. The signals are 

filtered and synchronously ,uadran»e demodulated, for example, usmg a Dolch- 
^"yshev m fiher. and dte complex signal values are d,en applied to cons»uct three 
matrices 410. 420, 430 which a« used to determine the mumal indactences and 
uormalizemeasurements.Ma,rix410isarawsignalmatrixmadeofthes,gna 

_ponen,sinaucedineachofa.ereceivercoiUbyeachofthe.ransm.nerco,,s.M^.^ 
42ol a 3X3 diagonal matrix, the three entries of which represent -he response of*e 
Ispec tive receiver coils to the caltWon coil. Matrix 430 is a 3X3 current matnx -th 
:resrepresent.g.ecurren.througheachco«(on.hediagonal)«.the.^^ 
„ 4,0 other coils (offfte diagonal). That is, it is a cross-couplmg matnx. TTeDSPor 
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processor also computes the inverse 435 of this matrix, and the current through the 
calibration coil is provided as a separate parameter 440. 

In general, the mutual impedance measuring system is designed for measurement 
of small mutual impedances at frequencies between 1 OOHz and 40kHz, and preferably 
between about 10 kHz and 15 kHz. The reference impedance Z.ef .may be a 10 ohm 
ilOOppm four-terminal resistor, so that is real and independent of frequency at the 
frequencies of interest Ihe whole system is ratiometric: the mutual impedances are 
measured with respect to Z^so that no precision voltage or current references are 
required. Moreover. Ihe frequency responses of the receiver measuring chamiel and the 
current measuring chamiel are Ihe same, so that factors of the complex gains that vary 
with frequency will cancel in the ratioing process. Similarly, frequency-independent 
complex gain factors of the receiver measuring channels v^ll also cancel in the 
determining the ratio of mutual impedances to T^. 

Overall, the system operates by converting the raw signal measurements as 
discussed above to a mutual inductance matrix, and then solving for position and 
orientation (P&O) from the mutual inductance. 

The first step to computing P&O is converting sum:of-products (SOP) data into a 
normalized raw signal matrix (RSM) of mutual inductances. Hie SOP data is processed 
to preserve phase information, so there is a real and an imaginary component of each 
parameter. TTie umiormalized matrix has these complex entries., while . the mutual 
inductance matrix consists solely of real numbers. As discussed above, it is Ihe end result 
of nonnalizing a complex signal matrix (CSM) for the effects of current, coil and 
calibration coil response. The goal of this step is to end up with a 3x3 matrix containing 
the mutual inductances between all possible combinations of sensor and transmitter coils. 
Because there are three sensor and three transmitter coils, there are a total of nine 
elements in the signal matrix. Each row of the matrix corresponds to a sensor coil and 
each column to a transmitter coil. 

The RSM is computed by performing a series of matrix multiplication operations, 
nae effects of transmitter current are normalized out (because the B field is directly 
proportional to the current through the transmitter), as are the effects of frequency (db/dt) 
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and receiver 
shown below 



.h^..gain(us^g*cca,«,«tioncon.«ponsen»*=s).T^f--ia» 



RSM ( Ae mutual inductance matrix) equals 



ICM 



(Fcal/FO 0 0 ) 
( 0 Fcal/Fl 0 ) 

(0 0 Fcal/F2) 



*CSM 

(Fcal_sop_C/Fcal_sop_0 ^ J^^ I ) 

*( 0 Fcal_sop_v. ^_ 1^^^ sopjC/Fcal_sop_2) 

( 0 

(cc_cflO] 0 0 ) 
* ( 0 cc_cf[l] 0 ) 
(0 0 cc_cf[2]) 

. ™ riF2 are Ae frequencies of ftex-.y-»d 3- transmitter coils, and 

whereF0,Fl,Rare4 ^ C is the measured sum^t-produc« cmren, fluough 

field, etc. 

disturbances). The processor ti • ^e done by pre-multiplying the RSM by 

:r:r::r— Lva.<.o,,«s.ne.c.is^ 

Appendix A. 



Specifically, 
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,h„e R is BOW a 3x3 martc wh«. .ffec« of re«iver and ^ans^tter gain and non- 
:Xona,i.esW^^ovoa.Oneoa.no«».n,a..e^.e..(.e^^^^^^^^ 
„f ^V^) and ft. magni^dc of *. :c..^ coordinates bas^d on ft= s<n«« of the «>.al 
field firom the three transmitter coils. 
The pertinent equations (1-4) are: 

r=* yl6.0*k^ /Stotal 
x = r* /k*yl(5.0*Sx-Sy-!iz)/y^ 



(1) 
(2) 
(3) 
(4) 



yzzr* /k*^l(5X)*Sy-Sx-Sz)/iri 

Z = r* /k*yl(5.0*Sz-Sy-^/i^ 
„herotisasc,efi>c»rfta.isft.prod„c.of*.axeaandnumberof»n.of^iver 

j::L«e,coi.=ndp«n,ea«H.,of^spaoe...o.nsftesnn,ofs.^^^^ 
c,enKn«s,«a»d.S.iy.&areftesn„.ofsigna.n««xe,en«n.ss,u.«dfo,fteX.Y.Z 

transmitter coils, respectively. 

B«=a«s.offteh«nisph.rio..an,Hgai.yassocia«dwiftft.-*er«.ep™ 
„3y always assun.ei.is operating in fte^henrisphe^. and singly take dotprodncu 

:leen Lsor ou.pu. vec«>rs (*e rows in a «ansposed signa, n.,rix, u. determ^e *e 

signof y andz. 

However, fte sohrtion is unsUble if any of fte coordinate values is close te ^. 
(See AppeodtxA). TO corr«. for ftis^eprooessor preferably ma—llyrou^ 

refers U. a region in which fte nvnredca, coordinate values are subs«n«ally non- 
^h:^chd«.!s«b.eso,u..,offtecoordinateva,uescanbeob.inedbyu«,^g 

r*ov:n.eftodoU,gy.1.erote.ionisdetern,inedbycon.pu.ngfte,ua~^^^^ 

te n,ove fte initial position estimate te an "ideal" solution location where 
:w:;l!ftcx.yand.a.es).™se„suresfta.fteso,utio„isnun,erica„yaccurateand 
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^^T,c rn (A'i with the rotated R matrix. After 
.ui Th^ processor then recalculates equations (IH^)Wixnui 

stable. The processor tnen ^^ocessor de-rotates with the inverse 

. . .1. 'WntpH" nosition coordmates, the processor 

e^loyed to d^tennine to sign of to coordinates. 

. . , and . values axe ton used to confute 3:t3 n.«iocs T. T,^^ 
,,e.es;eI.Hcesr«r.a,e.tationn,aMces^iU^to(p.sition.,)sp.^ca, 

coordinates ang.es. a. A -pecHvely. Note, however, totbeeanse 

cos(a)= x/sqrt(x^ + y^) 
sin(fi)=z/r 

... - /T«-' may be constructed directly from x,y. and z. 

To; T/Sani f/ieir inverses (Ta . ip J may 

inverse-coupling matrix, Q (descnoea in u 
Basically, 

(5) ff'.Ta-"Tr'*S-'*Tfi*Ta 

Wi«.ff'andtosi^a,n.«ix.«.toprocessoreon,p«testorota.ionn».rix.A.tha. 
rlazero—n sensor outputintotooutputoftotruesens^. 

(6) ^^Cr"/*)*!?'!?' 

„,er.toe,en,entsofAcanb.descnW*yarotation,uaten.ion.,-=ft.^.^.«»>: 
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A[0J[1 J= 2*(qy*qx+ qs* qz), 
A[0J[2J = 2*(qz*qx-qs* qyh 
A[1][0J = 2*(qy*qx-qs* qzh 
A[1][1J= (qs' - qx^ + qy^ - 
A[l][2]= 2*(qz*qy+qs*qx)' 
A[2][0J = 2*(qz*qx+qs*qyh 

A[2J[1] = 2*(qz*qy-q5*qx)' 

A[2J[2J= (qs^-qx'+qy^^qz'>- 

^eelen,entsof.*maybeusedtofindarotationquaten.ion'.%usi^^^^^^^^ 

* A' A^nf^ev;iX>CTofB.liom,Closed.form solution Of absolute 
p.ocedureinAppend«A8of1hepap^f ^ , Note, 

orientationusmgunitquatermons,J.OptSoc. Amer ^^...^ 
.owever^athere^islheinverseof^erotationmatnxde^T^mlh^^ 

fo.goingp«>cedureprovides a very good ""'^Z^^^^ 
orientationofthe sensor. Torefineihe estimate ofposidon and onent^^^^^ 

:Cesafitting.suc.asaf^t-orderorleasts^e.^^^^^ 

..e:aseto..ode.st.at...^^^^^^^ 
inductances in the RSM. Hie idea is to compute a set 
to the estimated P&O. The basic equation is: 



(7) 



E+ S* {deUa_P_0_vectoT) = R 



where 

that is to be solved for. jft,rhv 
•n,e fo«rt, compo^n. of *= de«a quaternion. Ms, may be solved for by 

noimalizing it to unit length. 
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Note that because the matnx.i IS noisquoic, ^ ^ 

W. Brognan. Modem Control Theory. Prentice-Hall. 1991. pp 
equation we end up with is: 

Finally, sensor position is adjusted by: 

x=x+Ax; y=y+Ay; z=z+Az 
ana the sensor orientation (as a quate^ion). .onenUs adjusts 

computed delta quaternion, qdelta : 
qorient = qdeUa* qonent 

nas. *= processor wi* SOP data from 4,e mag„e«c .racking uniB 

^ refine positior. a.,1 orienUtion of *e movable coil assemb^.. 

Asr^rea above, ccr^in — P-"«'"*'»"'=^'^"'7^L 

is addressed in one or more of several ways. 

Oneof*es.«aysis«>moun,*ereceiver(sensor)of*esys.emonorc,oseu>a 

.ajor in^rferin, s«cn„e. sue. as (in a medio.. eon«x.) .be C arm of a fl_ 
asslbly The ,ransmi»er eoil assembly is *en .be .^vable assembly of *e «a*er and 
be generally posi^onedremo^ly ,»om*ein«rferings«emre. Tbe n^nra, 
ITl Iker If L invention ean be eonfigured « opera» effective^, wft small 
":rb:"arHedby.bana-be,d.oo,and*elike.*erebyre„d^^^^^^^ 
a^Lcb feasible. P.r*er. .euse of mun.. ind»e.«.eebe.„een «,e — and 
liver „ni. in *e elec«omagnetie «aeker of O-e invention aUows s^appmg *e 
r „a,i.yof*e— randreeeivern.i«.e.g..utilizing.nominal^^^^ 

and/or a nominal «nsmi«er oo.-. as a receiver. "Hns recproerry ean be 
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"^cularty advantageous because, in some erases, itmay be more convenient to map 
and/or model the field generated by one unit, e.g.. receiver, rather than the other, e.g.. 



transmitter. 



Another aspect of the invention dealing with the field disturbance problem is to 
mount a standardized conductor about such an interfering structure to operate as a shield 
(e.g.. with respect to disturbances originating within the shield) and a known disturber, as 
seen by a sensor positioned outside the shield. The standardized conductor, such as a 
metal can. flien has a fixed position relative to one of the coil assemblies, and may be 
effectively modeled. Thus, the eddy currents induced in the sheet metal cylinder by the 
magnetic field from the transmitter assembly, and the secondary B field formed by these 
induced currents, may be modeled such that the processor can apply determinable field 
contributions from the shield. Without such a shield, applicant has found that when a coil 
assembly is employed in proximity to certain interfering structures (such as the imagmg 
assembly of a fluoroscope). the disturbance caused by that structure may vary 
considerably for two different but apparently identical mxits. The shield moderates such 
differences; electromagnetic fields originating within the shield are largely null outside 
the shield. Moreover, little of the tracking signal origination outside the shield penetrates 
inside so that the level of eddy currents induced in the structural conductive elements 
inside the shield are subject only to quite low fields and their eddy currents and secondary 
field effects are greatly diminished. These combined effects allow a tracking system to 
operate effectively by applying standardized (modeled or calibrated) corrections for the 
shield alone, without regard to the specific structure(such a s a fluoroscope imaging 
assembly) that it isolates. 

Preferably, the shield is made of a single piece of a highly conducting material, 
e.g.. a metal such as silver, copper or aluminum, or a conductive composite. The shield 
surrounds the targeted component or piece of interfering equipment, such as the image 
intensifier assembly of a fluoroscope. as much as possible, consistent with system 
constraints (e.g.. X-ray transparency of the top end. mechanical protrusion limits for 
camera movement, etc.). Preferably, the shield thickness is sufficient, e.g., more than the 
skin depth of the metal or more than several skin depths, so that the magnetic field it 
• generates will not appreciably penetrate the shield or extend into the volume behind the 
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jy^ field distortions introduced by *e shield may be mapped and a traclcer- 

* ^ •« , preliminary set-up or calibration step, or the effect 
correction table may be generated ma preliminary sei p 

of the shield may be modeled. 

evHi fei image intenrifiers of differait models. 

Figure 6 in„s«a.es such a sys.m in an «»bodimen. 10 use i. sa«i^ 
^vig. Jin an ope.«ng room en—. As sho«n. *e system 10 mo*^es a 
ZLcope20.aw«ksu6on30havingo„e<.rm<«aisplays32 33.nda 

^luse«roe.er„serin«rface34,an.aplurali.,of»ao,™ge«^^ 

T3 Tl«displaysmayillu«a«in,agedviewsorsynfl.es,z.dv,.ws40 oftool 

T3. ineuispiij J „.,nuinnslratedasaC-atmfluoioscopem 

mjecu^es. tissue or *e like. The fluoroscope 20 « .llusttated 

wl^h an x-ray source 22 is mounted on a s»ucu.ral member or Cann ^^^^ 

r„y receiving ^ detecdnguni.. referredto herein as an imaging assembly 24. Ibe 

rres ab^u. a patien. for producing .»o dimensional projection images of *e^ 

::en..mdi.eren.an.e. ^rarroZT. 

:::Lcr::::ir=ro:.^ 

«c,L,ormo.e«>epatien.maybemovab.ypos,noned.ln..«^^^^^ 
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^«iHe.oo,,sucha.aca.he.r.f.«ib..endo««pc»ra«ar«cuU«dU»L 
^ fcettaetoTl is preferably a smalUocW Clemen. pos,Bon=d m or a. fte 
cases, uic tia r _ , T1 » Trimire 6 to track coordinates 

operative tip of .he tool as shown by catheter tracker Tl tn Ftgurc 6, to tr. 

of the tip within the body of the patient. 

AS will be understood by those having ordinary sldB in *e art. fh.oK,scopes 
.ypieally operate wis, thex-ray source 22positionedoppositcthe^ja or i^ge 

!Lng aslbly 24. While in son» systems, the X-ray source is fixed overh«4 and the 

' . ^ . Hkcussion below will be illustrated with 

caineraislocatedbelowapatientsupport,thediscussionDe 

. t^miral r-arm fluoroscope, in which the source ana 

regard to Ihe more complex case of a typical C arm iiuor v 

2^areco„neo.edbyas.ructura.men,ber.*eC-arn,26.««.al.owsn»veme..of.he 
i:,nacan.era.ssen.blyabout.hepatientsoitn,aybeposi.ioned»prod^x-ray 

views fiom different angles or perspectives, these devices, the .ntagmg beam 
r^lyd.verses..anan.le,there,ativelc^tionsandoHen.atio.ofthe««^ 
can^raLywithpositSondnetostrucmral flexing andn>echan,canooseness.and.he 

^ o7bo* «,e source and the cantera With respect to the parent and/or a too, wh^ 
::desi..dto..ackn.ay-.va.yinditferentsho.s.Thus.theuseofoneorn.re^^ 

l.en.sT3(trans..«ersorrec.ivers)o„the structure itself ntaybeu^^^^ 
.^se n^ventents in the surgical navigation systen.. Fur^er detatls of *,s complex 

imaging environment may be found in commonly owned United State Patent 
Zl^onsSerialNos. 09.«,.940a„d0,/360,e0S.bo.filedonApr^l2S 2^^d 

elofwhichis incorporatedhereinbyreference. Bachof*osepat«^t^PP c^ ^ 
discloses a tracking system wherein tracking of the fluoroscope as well as the ««1 and 
patient produce image sets and surgical navigation of «Aanced accuracy. 

m imaging beam ilhts^ted by B in Figore 6 diverges ftom th. source 22 in a 
generally ^uncatedconicalheam shape, and theC-arm 26 ismovablealongagene^ly 

Lua.ea.h.opositionthescurc.andcamera^rimaging^mdif..»t — 

This generally involves posidoning the assembly 24 as dose as pcsstb^ behmd the 
relevL tissue or operadng area of ,h.patient,wHletheC.armassemb^ts mo^^^^ 
rou^ly about a targeted imaging center to the desired viewing angle. Tlte Carm o o*er 
ZL Lcture 26 m.y move eccen^cally or transla^onally to allow better clearance of 
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coordinates. 

^ K ™ i„att«=kerofft«i»vention.ashieldmaybemodeledin«l» 
AS discussed above. » » i„i«a..y. 

either obviatmg or reducing me orientation components by 

::r .te^^^.— ^^^^^^ 

flu„„,sc«pe™.ybee^-ey ^^^^^^^^ 

position. Tins may produce a kerne. . „ a fest or even closed form calculation 

.3Ssivestoredfieldoa.ibra«on,*les.orn«^yreM^^-^^^^^^^^ 

rrmCr:;:e«e.ddis«^.io...-*-— 
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^ analytic solution for position and orientation coordinates 

This assumption permits an anaiyi However, in cases 

^ • „w"P&0"^ based on the measured sensor data, nowc 

„,^.he ac^a. fields de^te^ P ^^^^^^^,^^^.^ 

eo« sues axe _ .enn se„s« is used 

^«tog or magneto — ^ ^ ^ig^,^ An 

,„aly«c sohnon for P&O » ™ g ^^^^ ^ 

-rirrd..esof.e— 

How^er. in acoordance ^ anofter aspect of the present Utventiot. *an 
JTw c-ting a distortion map for every application, a system -ntplentent^ a 
proceedmg by creatmg a u , „f the sensor field (for example, 

P&O dete.n»«tio.basedo„anK,re accurate model oftheset^^J^^^ 

modeling U,e extrinsic distortion, that eUminates the need for a map by cnr,ng 
distortion problem at its source. 

. • a, models (e g finite element methods) can be mad. arbitrarily accurate. 

.^....rz::::-*--''-^-— ^^^^^^ 

-^TT ™^Ibe^in.itedby the ability to match the actual geometry 

Wsapproachhasbeen applied inafin.h.x aspect ofthe invention toa 
„ 'clclcingsystemwhere*esensorismo«ntedclosetoaIargepara^,nc 
flucroscoptc tr^^^^ng sy ^ ^ ^ 

dtSnent.namelytheC-armora k „„^a«a ™ fte c-arm image 
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(i) 



(ii) 



axis dipole located near a filamentary conducting ring. Tbe fields of bolh a dipole and 
filamentary ring can be computed analytically (see. for example. Smylhe. WR. Static and 
Dynamic Electricity. ncmisr>h.r. Publishing, New York. 1989. A composite model of the 
coupling matrix components can therefore be written as Ihe sum of analytic expressions. 
The following sections describe the model in detail. 

To embed the new field model in the tracker fit routine, one needs a way to 
compute fliecouplingmatrixbetween the transmitter and receiver in the presenceofthe 

distortion source. To accomplish this, applicant makes two assumptions: 

Assume that the ixansmitting and receiving sensor coils can be modeled as 
mastic dipole elements. This assumption allows for a conv^ient dot 
product relationship between the field vector and mduced voltages, as the 
coils are assumed to be point sources/receivers. 

Assume that the dominant distorter is the metal can around the hnage 
^^Ser, and that we can model that can as a -7^; 
assmnption allows us to use analytic ' expressions for ^^ T^^f^^^ 
radiat^by the ring, which are available m *el^t«"t"^«(Smy^e, supra 
or JacksoX JD Classical Electrodynamics. Wiley, New York, 1975). 

using these assumptions, the development proceeds as follows. First, use network theory 
to derive an expression for the coupling between the transmitter and receiver m the 
presence of the ring. This expression contains several mutual and self-inductance terms 
Aat we need to compute. Known expressions for flie coupling between two dipole 
elements are used together with the fields generated by a ring to compute the required 
inductances. Finally, the signal matrix elements can be computed as a fiinction of P&O. 

The voltage induced on a receive coil due to a current flowing on a transmit coil 
in the presence of a metal-ring can be computed by considering Ihe system as a three-port 
network, as schematically shown in Figure 7. 

. The impedance parameters describing flie coupling among the ports are given by: 
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(See, e.g.. Pozar. DM. Mcrowave £„gmeerf«g.Addison-Wesley. Reading, 

Massachusetts, 1990.) 

Cc»«ider*e case where w.wan.«.taowft= voltage mduc«.onpon2(*. 

receiver, when po« 1 (*e ^«^r, is driven wia, a fixed vo,.ge. po« 2 . loaded w,d> 
: Uanoe Of Z. and port 3 (d.e ring) is .oaded wi* a aho« cU^it so^e 

:gehL.weeanwri»*ei.d»eedv.i«.geonpori.,eWve,o*ecn™n.throughportI 

as: 

. (Eq.5) 

; L -.. ^ . 

Kwe assume *a,Z. is "large" relative to dKO*er.e>ms in *«scco»Jsetof 
rackets, and that the resistance of the ring is negligible. Equation 5 reduces U>: 

y, , (i<ot.,.Xj'^2») 

■ ^ ■ i,'«iri,mt for many ciSes, and a>e full expression in (5) can 
This approximation is sufBcient lot many t--^. 

be implemented where the additional accuracy is required. 

Tocalculatethemnmalinductancehetween^etransminingorreceivingcoflsand 

thering. we apply reciprocity and letdieringact as ^csou,^ for both cases. Byfi^ 
ir^jfieldradiatedby.currentfiowingon^ering. the voltage mducedatd- 

Ismitt^receivecoilisgivenbyasimpledotproduct Under the ma^eo-stat^c 
::^tion,ti.eB.fieldprod„cedbyauniformc„irent,nowingonafilamenta,ynng 

(Figure 8) is given in cylindrical coordinates by Eq. (7): 



(a - p) + 2 J 



Lei. and^are complete elliptic integralsofthe first and second^^^^ 
and the argument k for the integrals is given by 
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102288-15 

4ap-^a + pf +z^y . 

Equation 5 requires the computation of one self inductance term and three mutual 
iiiductance terms (additional terms mustbe calculated for Equation 6. which makes no 
approximations). 

The self-inductance term is for the ring, and is given by [7] 



where 

a is the radius of the loop, 
b is the radius of tiie wire, 
fi' is the permeability of the wire matenal, and 

in [10] using the dipole assumption, and can be calculated as 
Atv ' R 

where . . 

A\ = effective area of fee transmitting dipole, 

ji, = unit vector in the direction of ^i, 
A\ = magnitude of y4i, 
Ai = effective area of the receiving dipole, 
A2 = unit vector in the direction of -^2, 
A2 =magnitude ofAi- 

^ = the vector from the transmitter to the receiver 
A = the unit vector in the direction of R 
R = the magnitude of R . 

Given a position and orientation of the sensor relative to the transmitter, Ae range 
R can be calculated. The effective area terms are calculated based on the coil geometry 

and a calibration procedure [1 1]- 
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remaintag ^vo n,»ma. inductance terms, ftose fc.^»«n .h. ttansmitting coi' 
a.d,h.ring,.ndbe,weenihe,«eivi„gcoila»d.he™g.n,ayb.cateu.a.edin.s.n„lar 



manner. 



The position and orientation of the ring relative «> the sensor coils is <o 
^ and feed. TT^erefore. to con.pn« -he mutual inductance L. in E„.tio. 6 firs, 

confute Ae f3.0 of the sensor in the coordinate 6»ne of the ring shov™ n. Ftpue 8. 
Next compute the magnetic field per unit current produced by fte ring a. dte ^ 
Zi^r. fU, *e mutual inductance term is calculated for each sensor col WJ^ 
the dotproductbemeen the coil vector and the field, using *e point source property of 

the dipole approximation. 

Computing fte mu«»l inductance be,«ee. the transmitter and the ring requires 
additional step. Giv«. the P40 of the sensor in the transmitter coordinate system, 
„efirs,conver.s..ha.w.havetheP&Oof*e.r.nsmitterinthesensorcoord,na.e 

ftame. T„e procedure then proceeds as above. The mutual inductance L„ ts computed 
by calculating titeP&Oofthe transmitter intireringcoordinate system. Therelationshrp 

between the ring and ftc sensor is imown and fixed as noted above, so once the 
nansmit.erP&Oistaownrelative.od.sensor,conversion.ransformattcnsa.eobv,ous. 
Next compute *.emagneticfieldpcruni.c„.rentproducedby*eringa.*ese^or 
location. Fin.lly.themu«.alinductancetenniscalcula.^ for each sensor coUby^g 
fte dotproductbetweenthe con vec»r and th. field, using the point source proper^, of 
the dipole approximation. 

Wim an of the induc^mc. terms available. Equation 6 is then readily calculated. 

overall. *e approach of measuring muhtal inductances aHows the system to d»el^ 
robustmeasuremen^inwhichrelianceonreferencessuchas^orast^alratioorus. 

of.common preamplifier with high precision output digitization, have elmtmated 
numerous sources of variability while allowing effective use of small magnetic 
assemblies to achieve accuracy and resolution. 

Moreover, by introducing a shield or virtual distorter and fixing a receiver 
assembly with respect to dte shield. d,e effect of individual distortion environments are 
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substantially eliminated, and the shield itself may be effectively modeled, eliminating the 
cumbersome requirement of compiling a mapping or calibration table. The model may 
then produce a seed calculation, allowing computationally effective fitting processes to 
replace the cumbersome calculations of the prior art. In other words, a distorter that is 
shaped to optimally shield magnetic fields from existing environmental distorters can be 
utilized without any regard for a specific desired field shape. The shield's effect on the 
field can then be modeled, or alternatively, mapped. 

The invention being thus disclosed and illustrative embodiments depicted herein, 
fiirther variations and modifications of the invention will occur to those skilled m the art, 
and all such variations and modifications are considered to be within the scope of the 
invention, as defined by the claims appended hereto and equivalents thereof. 

What is claimed is: 
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